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A new position-space renormalization group approach is investigated for bond directed percola- 
tion in two dimensions. The threshold value for the bond occupation probabilities is found to be 
Pc = 0.6443. Correlation length exponents on time (parallel) and space (transverse) directions are 
found to be !/|| — 1.719 and = 1.076, respectively, which are in very good agreement with the 
best known series expansion results. 
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I. INTRODUCTION 

Directed percolation has been studied extensively [1- 
6] , since it plays very important role in a large variety of 
non-equilibrium systems with a single absorbing state . 
A wide array of dynamic processes as fluid flow through 
a porous medium in a external field |^ , forest fires PpO) 
or epidemic growth models [prl , reaction-diffusion sys- 
tems [0[_^4|, damage spreading ||l^, self-organized crit- 
icality ||l6|] , models of growing surfaces with roughening 
transition []l7|-^, etc. fall into the same universality 
class as directed percolation. There is no exact result 
for the critical exponents characterising the dynamical 
phase transition seperating the absorbing steady state 
from the active phase. The correlation length exponents 
in the longitudinal and transverse directions, defined via 
^11 b — Pcl""^" and \p — Pc\~^^ and the percolation 
threshold pc have been found to great accuracy by series 
expansion methods ||,^ to be v\\ = 1.733, v±^ = 1.097 
and Pc — 0.6447, respectively. 

The fixed-scale transformation (FST) introduced 
by Pietronero, Erzan and Evertsz ||2|,|2^ has been used 
to investigate the self-similar properties of irreversible 
growth models, namely, diffusion limited aggregation ||2^ ] 
and dielectric breakdown |24[ , as well as other intrinsi- 
cally critical growth problems such as cluster-cluster ag- 
gregation . The FST has also been applied to systems 
such as percolation [ 26 1 , Ising or Potts models |27j , inva- 
sion percolation ||2^,2£| and to directed percolation | |3C| ] 
at the critical point. Renormalization group ideas have 
been used in conjunction with the FST approach | |3l[ |, in 
order to identify the scale invariant growth rules in terms 
of which the finite cell fixed scale transformation matrix 
should be computed, with good results. Forest fires | |3^ ] 
(which fall into the same universality class as directed 
percolation) and the sandpile model |3^,^^ have been 
studied by means of the "dynamically driven renormal- 



ization group" pSq l , which combines real space renormal- 
ization group ideas with a dynamical steady state con- 
dition reminiscent of the fixed scale transformation ap- 
proach. 

In this paper, we would like to introduce a novel po- 
sition space renormalization group (PSRG) treatment of 
the directed percolation problem. This approach mod- 
ifies conventional PSRG methods in two ways: i) the 
weights of different initial states in the RG cell are com- 
puted from the steady state distribution found from the 
fixed point of the FST; ii) a. "dynamical" coarse graining 
procedure is defined which allows for the appearance of 
two different scale factors in the longitudinal (time-like) 
and transverse (space-like) directions, thus taking into 
account the self-affine nature of the problem. These scale 
factors are determined independently, without having to 
make any additional assumptions. 

The paper is organized as follows. In the next section, 
we construct the renormalization group transformation, 
subject to the dynamical fixed point condition, and com- 
pute the critical parameter as well as the transverse and 
longitudinal correlation length exponents. In the last sec- 
tion we provide a short discussion. 



II. RG TRANSFORMATION WITH DYNAMICAL 
COARSE GRAINING 

We consider 1-time and 1-space dimensional bond di- 
rected percolation as a growth process. Let us recall the 
growth rules: if a site at time t is active, its two neigh- 
bours at time t-\-\ may be activated, each independently, 
with a probability p. If p is larger than the threshold 
value Pcj there is a finite probability that the growth 
process is continued indefinitely. 

We now outline a renormalization group procedure 
which takes into account fluctuations in regions larger 
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than the renormalization group cell by using a steady- 
state distribution of initial conditions, obtained from 
the FST. It also introduces self-consistently determined 
rescaling lengths, rather than pre-set scale factors be- 
tween the original and coarse grained lattices. 

Our renormalization group cell is shown in Fig. 1(a). 
The boxes AA'BB' and CCA' A" will coarse grain, un- 
der a dynamic coarse graining procedure we describe be- 
low, to the bonds ab and ac, as shown in Fig. 1(b). This 
process conserves the transverse and longitudinal direc- 
tions, but does not conserve the lattice angles, since these 
two directions are scaled differently. Rescaling by the ap- 
propriate longitudinal and transverse scale factors will re- 
store the original lattice (Fig. 1(c)). The lattice spacing 
is taken to be \/2 for convenience, yielding £± = 1 and 
£\\ = 1 for the transverse and perpendicular distances 
between the nearest neighbors, such as AA' and BB' . 
Clearly, we only have to consider one of the boxes, e.g., 
AA' BB' for our renormalization procedure, by symme- 
try. 
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tial configurations have to be considered as illustrated in 
Fig. 2. Any spanning configuration may have more than 
one path which connects the origins to the end points. 
For instance, two different paths are possible for con- 
necting the sites A and B' in the spanning configuration 
shown in Fig. 3. 
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FIG. 2. Four different initial configurations are repre- 
sented. In computing their respective weights, (see text) the 
bold lines in the each configuration are replaced with the prob- 
ability p, and the dashed lines with 1 — p. 

In order to compute the probabilities of finding each 
initial configuration, say Wi, i = 1, 4, we make use of 
the time invariant probabilities for the relative frequency 
of doubly or singly occupied boxes in a box covering of 
the equal-time transverse subsets on the cluster of active 
sites, in terms of the probabilities p. 
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FIG. 1. The dynamical coarse-graining procedure, (a) The 
RG cells AA'BB' and CC'A'A" (the bold lines) in the original 
lattice are coarse-grained (fe) to the bonds ab and ac respec- 
tively, which are rescaled in the next step (c) to preserve the 
lattice angles. 

To obtain the renormalization transformation for the 
bond occupation probability p, we compute the total 
probability P{p) that a spanning path exists across the 
box AA'BB' , by considering all the different initial con- 
figurations (states of A and A') and the spanning con- 
figurations that can be obtained from them. A spanning 
configuration is defined as one where a path starts from 
either A, or A' or both, and ends on B or B' or both. 



A. Steady state distribution of initial conditions 

Since there are two possible origins which can be ac- 
tive both together or by themselves, four different ini- 




FIG. 3. Two different paths belonging to the same span- 
ning configuration. 



In recent work |30|, the FST approach has been used 
to calculate the fractal dimension of directed percolation 
at the critical point, given the threshold value of p. The 
main idea underlying the FST approach here is that at 
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the critical point, as t —> oo, the infinite transverse sub- 
sets of active sites at any given t are statistically similiar 
under translation in t. In particular, they can be modeled 
by Generalized Cantor Sets generated by a random se- 
quence of fragmentations obeying a one-parameter, scale 
invariant distribution as shown in Fig. 4. Given a hierar- 
chical partitioning of the transverse subsets into cells of 
size 2~^, the relative probability of encountering singly 
occupied or doubly occupied cell among nonempty cells 
at a fixed arbirary scale is equal to C or 1 — C. This 
probability can be computed |Q from the fixed point of 
the fixed-scale transformation as. 



C 



2p2 + 4p _ 4 



(1) 



up to first order in the FST approach. 
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B. The renormalization transformation for p 

Depending on the initial configurations i = 1, 4, the 
total number of possible spanning configurations will be 
different. One sees that in Fig. 2 (1) only one; (2) seven; 
and (3) — (4) eighteen spanning configurations each are 
possible. It should also be noticed that some spanning 
configurations can be observed in more than one different 
initial configuration. For example, the spanning config- 
uration which contains only one path which starts from 
A' and ends at B' can be observed in all the initial con- 
figurations except the second one (see Fig. 2). The total 
probability fi of the spanning cluster for the i'th initial 
configuration is given by 



f2ip)^P^+p'~P^ 
/3(p) = 2p2+p3_3p4 

h{p)^h{p) ■ 



(6) 
(7) 
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The renormalization group transformation for p is then 
found to be 

P{p) = Wihip) + W2h{p) + W^hip) + Wifi{p) = p' . 

(10) 

The fixed point of this transformation gives the thresh- 
old value, Pc — 0.6443 which is in agreement up to the 
third digit with the series expansion 1^,^ result, namely 
0.6447. 



C. The Affine transformation in the longitudinal and 
transverse directions 



FIG. 4. The generator for the random Cantor set corre- 
sponding to a transverse subset of the incipient infinite clus- 
ter. Cells of type-1 and type-2 participate in the fine graining 
process with corresponding probabilities Ci and C2 . We rep- 
resent the activated sites with black circle and inactive sites 
with white circle. 

We now make use of this FST results to compute the 
normalized weights Wi of the initial configurations i in 
terms of C and p, 

Wi = pC/2 (2) 

W2 = {l^p)C /2 + {\-pf{\-C) (3) 

W^3=p(l-C/2) (4) 

Wi^ p{l^p){l-C) . (5) 

In the following subsection we proceed to obtain the 
renormalization transformation in terms of which the 
critical value of p can be determined. 



The system has two independent correlation lengths ^|| 
and ^j^, parallel and perpendicular to the time direction 
respectively, which diverge with different exponents as 
C|| ~ |p — Pcl^''" and ~ \p — pc\~''^ ■ To compute these 
exponents, we use the eigenvalue equations. 



dp 



= b 



II 



(11) 



dp' ^ 

p—Pc " dp p—Pc ^ 

We determine the appropriate rescaling lengths 6|| and 
b± from, 



^11 = 



< i|| > 



In 
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in terms of the average projected lengths, < L|| > and 
< L± >, of the spanning paths onto the time and 
transverse directions as shown in Fig. 5. These quanti- 
ties are the amounts by which the coarse grained lat- 
tice (see Fig. 1(b)) has been dilated in the "dynamical 
coarse-graining" step. Since we have taken the lattice 
constant to be v^, the projections of a single bond on 
the original lattice onto the time and space directions are 

e± = h = 1. 
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III. DISCUSSION 






FIG. 5. All possible paths and their extremal projections 
contributing to the dynamical rescaling factors < > and 
< L_L >. See Fig.l. 

In the time direction, one has to take into account 
the fact that under coarse graining, different time steps 
collapse onto each other; that is why we consider each 
path originating from A or A' and terminating on B or 
B' as contributing seperately to While enumerat- 
ing the possible paths over which the average is taken, 
it makes a difference if there is a bond between A and 
A' , since the paths change in case A' is activated by A. 
The same thing holds for the end sites B and B' . Note 
that the result of taking the extremal projections first 
and then averaging is different from finding the "average 
path" and then taking its projections. This point will be 
further discussed below. 

For X = {_L, II}, we have 



There have been some earlier studies of directed (there- 
fore self-affine) systems via position space renormaliza- 
tion group (PSRG) techniques. Introduction of two dif- 
ferent scaling factors 6|| and b± due to existence of two 
independent scaling directions was firstly suggested by 
Dhar and Phani |^^, where they employed a decima- 
tion transformation. However, their results are far from 
the accepted values of the critical exponents and depen- 
dent upon their choice of the RG cell, which determines 
b\\/b±. Very large-cell PSRG calculations by Zhang and 
Yang 1^ are able to accurately reproduce the self-affine 
behaviour of directed self-avoiding walks, but are more 
appropiately in a class with Monte Carlo RG. A bond- 
moving and decimation transformation for anisotropic di- 
rected bond percolation in arbitrary dimensionj 38| and 
its generalization for other directed systems [^| by da 
Silva and Droz give the critical fugacity very accurately 
in all dimensions for directed self-avoiding-walks. For 
two dimensional directed percolation, it yields good re- 
sults for the threshold value Pc and , but is not as good 
for v± . The approach in these papers is in fact very close 
to ours; however, the way the projections of paths are 
taken to compute the scaling factors in the two different 
directions are different. 

One may consider two different ways of projecting a 
path. The approach of da Silva and Droz [^Sjjs^] is to 
draw a vector between the origin and the end point, and 
take its two components to be the projections of the path. 
Our is to take the projection to be the transverse or longi- 
tudinal distance between the extremal points of the path, 
as shown in Fig. 5. 



X 



(13) 



/3 



where qi^p is the relative probability of any spanning path 
/3 belonging to an initial state i, and is to be found from 
fi by giving equal weights to the distinct paths in any 
spanning configuration equally. Substituting the value of 
Pc found from Eq.(IO) into Eqs.(2,3), we find. 



and 



< L_L 1.7996 



< L|| >= 2.5561 



(14) 



(15) 




These values yield, together with Eq.(ll), the correla- 
tion length exponents to be ly^ = 1.719 and vj^ — 1.076 
which are comparable with the best known results 
I'll = 1.733 and = 1.097. 



FIG. 6. Macroscopic paths with identical end-to-end vec- 
tors, but strongly differing extremal projections. 

We believe that our method yields a result which cor- 
responds more closely to what is meant by the size of the 
transverse and longitudinal fluctuations, in that it mea- 
sures more accurately the actual size of the region over 
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which a coherent flow takes place. This becomes more 
evident if we consider the incipient infinite cluster of our 
system for p ^ pc- In this cluster, there can be many 
paths which connect some point at time to to an end 
point at to + t- Even though the total number of bonds 
in these different paths are the same, the spatial size of 
the spanned regions can be radically different (see Fig. 
6). Since the transverse (longitudinal) correlation length 
corresponds to the spatial (time-like) size of the fluctua- 
tions, it is more appropriate to take the average over the 
extremal transverse (longitudinal) extent of each path. 

In conclusion, by incorporating the FST fixed point 
condition in the determination of our distribution of ini- 
tial configurations, and by a new, dynamical coarse grain- 
ing procedure which makes use of averages over extremal 
projections of spanning paths, we have succeeded in com- 
puting the percolation threshold and the correlation func- 
tion exponents much more accurately than before. 

Note added in proof: D. Sornette has kindly brought 
to our attention Refs. where the equivalence of 

directed percolation to a parallel version of the Bak- 
Sneppen model | p2| has been rigorously demonstrated. 
Also see for a comparison of directed percolation with 
sandpile models. 
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